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Optical properties of semiconductors can exhibit strong polarization dependence due to crystalline 
anisotropy. A number of recent experiments have shown that the photoluminescence intensity in free 
standing nanowires is polarization dependent. One contribution to this effect is the anisotropy of the 
dielectric function due to the fact that most nanowires crystalize in the wurtzite form. While little 
is known experimentally about the band structures wurtzite phase III-V semiconductors, we have 
previously predicted the bulk band structure of nine III-V semiconductors in wurtzite phase. Here, 
we predict the frequency dependent dielectric functions for nine non-Nitride wurtzite phase III-V 
semiconductors (A1P, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs and InSb). Their complex dielectric 
functions are calculated in the dipole approximation by evaluating the momentum matrix elements 
on a dense grid of special k-points using empirical pseudopotential wave functions. Corrections to 
the momentum matrix elements accounting for the missing core states are made using a scaling 
factor which is determined by using the optical sum rules on the calculated dielectric functions for 
the zincblende polytypes. The dielectric function is calculated for polarizations perpendicular and 
parallel to the c-axis of the crystal. 

PACS numbers: 71.15.Dx, 71.20.-b, 71.00.00 



I. INTRODUCTION 

Frequency-dependent reflectivity and absorption spec- 
tra have been a fundamental tool for studying the band 
structures of semiconductors. In semiconductors with 
wurtzite (WZ) structure, the optical properties are very 
strongly polarization dependenP"^ due to the crystalline 
anisotropy. The optical selection rules for WZ^l dic- 
tate which polarization dependent interband transitions 
are allowed. This helps in further resolving the band 
structure for WZ and is used to determine crystal-field 
splitting and spin-orbit splitting energies. 

A number recent experiments have measured the op- 
tical polarization dependent photoluminescence (PL) in- 
tensity in free standing nanowires^M^. These results re- 
veal a very strong polarization dependence in the PL in- 
tensity depending on whether the incident light is po- 
larized parallel and perpendicular to the growth axis of 
the nanowire. This polarization anisotropy has been at- 
tributed to effects such as heavy and light hole mixing^, 
large dielectric contrast between the nanowires and its 
surroundings^ and intrinsic band structure properties of 
the nanowire^. 

However, it is well known that the nanowires tend 
to crystalize in WZ phase rather than zincblende(ZB) 
(which is the more stable phase for non-Nitride bulk III- 
Vs). The polarization dependent optical anisotropy of 
the nanowires will have a contribution from the WZ crys- 
tal structure. However, without bulk samples it is not 
possible to determine the optical properties with tradi- 
tional reflectivity measurements. 

In Ref. [T5] we predicted the bulk band structure 
of nine non-Nitride III-V semiconductors in WZ phase 
using empirical pseudopotentials, including SO interac- 
tions. These calculations were based on transferable 



model pseudopotentials assuming ideal relations between 
the ZB and WZ phases and their lattice constants. The 
spherically symmetric ionic model potentials for the ZB 
phase were first obtained by fitting the calculated transi- 
tion energies to experimental transition energies at high 
symmetry points. The WZ phase band structures were 
then obtained by transferring the model ionic pseudopo- 
tentials to the WZ pseudopotential Hamiltonian. This 
is justified because in both ZB and WZ, all of the near- 
est neighbors and nine out of the twelve second nearest 
neighbors are at identical crystallographic locations^, 
while all the second nearest neighbors are equidistant. 
This method has proven to be quite successful in the 
past, in obtaining band structures of polytypes^HH] Our 
agreement with experiment was excellent for the known 
band gaps of InP, InAs and GaAs. 

In this article, we predict the frequency dependent lin- 
ear dielectric function for the WZ phase of nine III-V 
semiconductors, based on our previous empirical pseu- 
dopotential band structure calculations. The dielectric 
functions are calculated in the linear optical response 
regime within the electric dipole approximation. The 
required momentum matrix elements are obtained by 
using the calculated wave functions from our empirical 
pseudopotentials^. Momentum matrix elements must 
be corrected because the pseudo-wave functions do not 
include the core states. Such corrections have been made 
using nonlocal correction terms^^l. Since the core 
states of the two polytypes should be very nearly the 
same, such core corrections should be the same for ZB 
and WZ forms. Therefore, we use the optical sum rules 
to obtain a set of normalization constants which yield 
eo = e(w = 0) for the ZB phases that fit experimental 
values. These normalization constants are then trans- 
ferred to their respective WZ polytypes to obtain their 
dielectric functions. 
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Our paper is organized as follows. In Sec. [TTJ a brief 
theoretical background of the optical properties is pre- 
sented along with our method for the calculations. The 
calculated dielectric functions and reflectivity spectra are 
presented in Sec. |III| along with our tabulated results for 
eo- This is followed by a brief discussion and summary 



II. OPTICAL PROPERTIES 

Consider a semi-infinite crystal having symmetry 
equivalent to or higher than that of the orthorhombic 
crystal system. If we choose our coordinate such that 
the z-axis is the surface normal and the x-z plane is the 
plane of incidence, then the reflectivity for light po larize d 
perpendicular to the plane of incidence is given b y 29 l 3 ° l 
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where at an angle of incidence. Similarly the reflectivity 
for light polarized parallel to the plane of incidence is 
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where rij (i = x,y or z) is the complex index of refraction. 
In case of an optically uniaxial crystal such as WZ, if the 
surface normal is parallel to the c-axis then n z = nn and 
n x = n y = n -L- In the linear response regime, n(uj) — 
\J e(w), which can be written as 



e(w) = e» + ie"(w). 



(3) 



In the linear response regime, the real and imaginary 
parts of the dielectric function are not independent of 
each other, but obey the Kramer-Kronig (KK) relations 
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where V is the Cauchy principle value. 

To calculate the optical properties of the WZ semi- 
conductors, we first evaluate e"(ui) based on our earlier 
empirical pseudopotential band structure calculations^ 
and then obtain e'(w) using the KK relations. In the 
electric dipole approximation, assuming only direct band 
to band transitions are allowed between an initial state i 
and a final state, j, e"(w) is given by 



7T e 

t 2/, .2 



(6) 



/ iMijfSiE^k) - E Vji (k) -?iLj)d 3 k 



where J BZ is an integration over the entire Brillouin zone 
(BZ), Y^jj is a sum over all initial valance band and final 
conduction band states, and E v (k) and E c (k) are the 
valance and conduction band energies at their respective 
ks. For the delta function, we use 



S(AE-huj) 
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where 7 is an adjustable damping parameter that can be 
used to incorporate lifetime broadening effects. We used 
7 = 100 cV _1 which gives a transition linewidth of about 
35 meV^. 

The momentum matrix elements Mjj for intcrband 
transitions are obtained from the pseudo-wave functions 
from our empirical pseudopotentials for WZ 1 ^, given by 



c (r)=^c(k,G)exp[j(k + G).r] 



(8) 
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where c(k, G) are the eigenvector coefficients at a given 

k. The momentum matrix element between the states i 
and j is 

M«(k) = mm) (9) 

where p is the momentum operator. Using Eq. (JsJ , 
Mij can be rewritten in terms of the expansion coeffi- 
cients as, 



My(k) 



G 



c*(k,G)c,(k,G)[(k + G)-e] (10) 



where e is the polarization vector. The expansion coeffi- 
cients, Ci and Cj, are the eigenvectors of i th and j th state 
in the pseudopotential Hamiltonian. e"(w) is then calcu- 
lated using Eqs. [5J [7j [10] for light polarized parallel and 
perpendicular to the c-axis. In practice, it is difficult to 
explicitly evaluate the Brillouin zone integral in Eq[7] be- 
cause of the prohibitively large number of k values that 
would be required. However, integration schemes that 
allow the BZ integral to be replaced by a sum over a 
set of special k points can be used. We have used a 
set of 4.5 x 10 4 special k points based on the scheme of 
Monkhorst and Pacl<P. 

Momentum matrix elements calculated using the 
pseudo-wave functions must be corrected since the 
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FIG. 1: Comparison between calculated and measured 
reflectivity spectra at normal incidence for zincblende 
phase (a) GaAs and (b) InP 
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pseudo-wave functions do not include the core states. 
One method is to include the commutator of the nonlo- 
cal pseudopotential and the position operatoi^lMI ) w hile 
Kageshima and Shiraishi have proposed correcting the 
momentum matrix elements by including a core repair 
terrrPsl. Both techniques cause small changes the dielec- 
tric function. Monachesi et al. have compared the di- 
electric functions of GaAs calculated with ab initio pseu- 
dopotential wave functions against calculations with true 
electron wave functions and find virtually no differences 
due to the missing core states^. In our calculations we 
take advantage of the fact that our pseudopotentials are 
being transferred between polytypes, and hence any core 
corrections should be nearly identical. We normalize the 
calculated e(uj = 0) to the experimentally known static 
dielectric constant by making use of the optical sum rule 



e = 1 + C- 



e"{uj)doj 



(11) 



where C is a scaling constant which is determined as fol- 
lows. First the dielectric functions for all the ZB phase 
semiconductors are evaluated using empirical pseudopo- 
tential wave functions. The constant C is then adjusted 
so that the calculated eo for ZB matches experimental 
values obtained from Ref. 1351 By empirically fitting to 
eo this method results in good agreement between theory 
and experiment for frequencies uj > 0. As an example, we 
show a comparison between the calculated and measured 
reflectivity spectra at normal incidence for cubic GaAs 
and InP in Fig. [I] The measured dielectric functions for 
GaAs and InP were obtained from Ref. [31)1 

For the WZ calculations, e" and ejj obtained using 
Eq. [7] are multiplied by the same scaling constant, C. 
We expect this method to yield good results for the WZ 
semiconductors since both polytypes consist of the same 
atomic species, so the corrections to account for missing 
core states will be the same. 



III. RESULTS AND DISCUSSION 
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FIG. 2: Real and imaginary parts of the complex 
dielectric function as a function of incident photon 
energies and polarizations for (a)E± in AIP, (b)i?|| in 
AIP, (c)E_l in AlAs, (d) E\\ in AlAs, (e)E 1 _ in AlSb 
and (f)£|| in AlSb 



Material 


eo 1 


e o 


AIP 


10.464 


8.232 


AlAs 


10.853 


9.165 


AlSb 


12.056 


11.933 


GaP 


11.708 


10.223 


GaAs 


12.481 


13.066 


GaSb 


15.215 


17.621 


InP 


12.812 


10.435 


InAs 


16.782 


13.610 


InSb 


16.952 


19.379 



The calculated eo for light polarized parallel and per- 
pendicular to the c-axis are listed in table jl} The values 
listed in the table show that generally, materials with 
heavier elements have larger static dielectric constants. 
We also see that ej~ < eJJ in the case of GaAs, GaSb and 

InSb, whereas e^ > e ' for AIP, AlAs, GaP InP, and InAs. 

For WZ materials, with the inclusion of spin-orbit in- 
teractions, all zone center states belong to either F 7 , r 8 
or Tg representations (all in double group notation) . The 
zone center conduction band minima have either r 7 or 
Tg symmetry and the valance band states are the T 9 
heavy-hole, r 7 light-hole and a Tj split-off-hole. The 
interband transitions for WZ are dictated by the polar- 
ization dependent optical selection rules^H For light po- 
larized parallel to the c-axis (-E^i), only transitions be- 
tween states with the same symmetry are allowed, i.e. 



TABLE I: Calculated static dielectric constant for light 
polarized parallel and perpendicular to the c-axis, for 
nine WZ phase III-V semiconductors. 



r 7 o T 7 , T 8 «-» r 8 and r 9 «r 9 . For light polarized per- 
pendicular to the c-axis (E±), the allowed transitions are 

r 7 o r 7 , r 8 o r 8 , r 9 o r 7 and r 9 o r 8 . Note, that 

the r 7 O r 8 transition is forbidden for all polarizations. 

The other high symmetry k-point transitions allowed 
for E\\ are A 7iS o A 7<8 , A g o Ag, K 4 r o #4,5, 
K 6 o K 6 , i7 4i5 -h- i/ 4 , 5 and Ha o H 6 . For E± the 
allowed dipole transitions are A 7 , 8 o A 7 ^, A 9 f-> A 7 , 8 , 
Ki^ o K 6 and H ii5 «-> H 6 . All M and U valley transi- 
tions (i.e. M 5 o M 5 and L 5 O L 5 ) are allowed for all 
polarizations. For a comprehensive list of optical selec- 
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FIG. 3: Real and imaginary parts of the complex 
dielectric function as a function of incident photon 
energies and polarizations for (a)E± in GaP, (b)^] in 
GaP, (c)£l in GaAs, (d) £j| in GaAs, (e)^ in GaSb 
and (f)£|| in GaSb 



tion rules for WZ at various high symmetry points and 
along various directions, see Refs. I5l and l37l . 

The zone center transition energies for the first thirteen 
states and their respective irreducible representations are 
listed in ReffTSl Our band structure calculations show 
that the indirect gap ZB semiconductors (A1P, AlAs, 
AlSb and GaP) have direct band gaps in WZ phased 
with Tg conduction band minima. These materials will 
be partially optically active, with transitions from the 
Tg heavy- hole bands allowed for E±, making them po- 
tentially technologically important. The other semicon- 
ductors with Tg conduction band minima in WZ phase 
are GaAs and GaSb. The Indium containing WZ semi- 
conductors all have direct band gaps with Tj conduction 
band minima and are optically active for all polarizations 
of the electric field. 

The calculated e'(uj) and e"{uS) for light polarized par- 
allel and perpendicular to the c-axis for the nine III-V 
semiconductors of interest are shown in Figs. 2]|4 The 
corresponding reflectivity spectra for both polarizations 
is shown in Figs. [5]|4] The reflectivity spectra for the dif- 
ferent WZ crystals are distinct and depend on the details 
of its electronic band structure. For the WZ semiconduc- 
tors, the most prominent features are typically seen up 
to about 7 eV. As expected, all the WZ phase semicon- 
ductors exhibit optical anisotropy, the degree of which 
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FIG. 4: Real and imaginary parts of the complex FIG. 6: Calculated reflectivity spectra at normal 

dielectric function as a function of incident photon incidence for (a) GaP, (b) GaAs and (c) GaSb 

energies and polarizations for (a)^^ in InP, (b)^ in 
InP, (c)£j_ in InAs, (d) E\\ in InAs, (e)E± in InSb and 
(f)£|l in InSb 



t; 0.6 
1 0.4 

I 02 





R i 




R || 


AIP 





5 10 15 

hw (eV) 



20 



0.6 



> 



8 0-4 
I 0.2 





R i 








R ll 


AlAs 





5 10 15 

hw (eV) 



20 



g 0-4 







AlSb 





5 10 15 

hw (eV) 



20 



FIG. 5: Calculated reflectivity spectra at normal 
incidence for (a) AIP, (b) AlAs and (c) AlSb 
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FIG. 7: Calculated reflectivity spectra at normal 
incidence for (a) InP, (b) InAs and (c) InSb 
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varies with the semiconductor. The reflectivity spectra 
for the two different polarizations shows several peaks 
which originate from interband transitions along various 
high symmetry points. 

In order to illustrate the the dielectric function's vari- 
ations about the fundamental absorption edge (FAE), we 
fit the numerically calculated e± and en (in the vicinity 
of their respective band gaps) to the Lorentz oscillator 
model. The real and imaginary parts of which are 

'M-Wg ^.q),.'^), (12) 

where, /, Qj and Tj were used as fitting parameters and 
are listed in table [IT] Note that the real and imaginary 
parts of the dielectric functions were separately fit to 
Eq[l2] and Eq{l3} 

It is seen, that in general, e" is greater than ej|. Based 
on optical selection rules, this can be explained by the 
fact that for E± a transition between the Tg heavy-hole 
(HH) and the T7 (or r 8 ) conduction band (CB) is al- 
lowed. Hence the FAE is at the bandgap (E g ). For E\\ 
the allowed transition with the lowest energy will be be- 
tween the F7 light-hole (LH) and T 7 . As a result for E\\, 
InP, In As and InSb, the FAE will be at Eg + Ax where Ai 
is the energy difference between the Tg heavy-hole and 
T-? light-hole. In general, Indium containing compounds 
have prominent features in the reflectivity spectra near 
the FAE as a number of valance band to conduction band 
transitions are allowed. In all others (A1P, AlAs, AlSb, 
GaP, GaAs and GaSb) no noticeable sharp features are 
seen in the vicinity of the FAE since they all have T s CBs 



and the optical selection rules forbid transitions from the 
LH and the split-off hole for all polarizations and from 
the HH for E ]v 

The spin-orbit coupling can alter the ordering of the 
valence band states in different WZ semiconductors. In 
our calculations, the valance band states in all materi- 
als except InSb have normal ordering (LeTg, Tr, T ) 7 ^ 38 ^). 
The ordering of the valence band states in InSb is compli- 
cated by the very large spin-orbit splitting which forces 
the F7 split-off hole bellow the next Tg state, resulting 
in the unusual Tg, T?, Tg, T7 valence band ordering. At- 
tempts to optically measure the spin-orbit splitting en- 
ergy in WZ type InSb could be complicated by this. 



IV. SUMMARY 

In summary we have calculated the optical proper- 
ties of A1P, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs 
and InSb in WZ phase using empirical pseudopoten- 
tials. Their complex dielectric function was evaluated 
as a function of incident photon energy up to 20 eV for 
light polarized parallel and perpendicular to the c-axis. 
The dielectric functions were calculated in the dipole ap- 
proximation by evaluating the optical momentum ma- 
trix elements using empirical pseudopotential wave func- 
tions (from Ref. 18) on a grid of about 4.5 x 10 4 special 
k-points. Corrections to the pseudo-momentum-matrix 
elements to account for the missing core states are in- 
troduced via a scaling factor, which is determined from 
the ratio of the calculated to measured static dielectric 
function for the corresponding zincblende polytypes. The 
reflectivity spectra for all nine WZ semiconductors was 
also calculated for both polarizations and was seen to ex- 
hibit optical anisotropy as well. We have also calculated 
the static frequency dielectric constants. 
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